Many-polaron problem by cluster perturbation theory 
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Abstract 

The carrier-density dependence of the photoemission spectrum of the Holstein many-polaron model is studied using 
cluster perturbation theory combined with an improved cluster diagonalization by Chebychev expansion. 



Key words: elcctron-phonon interaction, cluster perturbation theory 
PACS: 71.27. +a, 63.20.Kr, 71.10.-w, 71.38.-k 



Polaronic quasiparticles originating from strong 
electron-phonon interaction play an important role 
in polar solids, including alkali halides, transition 
metal oxides/perovskites and, in particular, quasi one- 
dimensional (Id) materials such as MX chains, conju- 
gated polymers or organic charge-transfer complexes. 
However, the theoretical description of the underlying 
models represents a challenging open problem. 

Recently, significant effects of finite carrier density 
on the spectral properties of a polaronic system have 
been discovered for the spinless Holstein model [1], re- 
vealing the shortcomings of widely used single-polaron 
theories. The Hamiltonian of the Id tight-binding (t) 
Holstein model reads 

H = c\ Cj + w b\h - /Bp^X; ^ + b i ) 

(M) » i 

Here cj (&|) creates a spinless fermion (a phonon) at 
lattice site i, and hi = c\c i . The parameters are the 
adiabaticity ratio UJo/t, and the dimensionless coupling 
constant A = E p /2t, with the polaron binding energy 
E p . 

While existing work based on finite-cluster calcula- 
tions was restricted either in energy or momentum res- 
olution [1,2], the cluster perturbation theory (CPT, see 
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[3] and references therein) yields accurate results in the 
thermodynamic limit. 

Here, we calculate the zero-temperature Green 
function G(k,uj) = {(c fc ; c|.))a, using CPT in combina- 
tion with a Chebychev expansion technique [4]. The 
latter does not suffer from the shortcomings of other 
methods, and exploits the analytical separation of the 
symmetric phonon mode [5] in order to ensure well- 
converged results. Calculations have been done on 
parallel supercomputers, with the largest matrix di- 
mension exceeding 10 10 . Within CPT, the one-particle 
Green function is obtained from G _1 = G' c ^ _1 — V, 
where G' c ' and V denote the corresponding real-space 
cluster Green function and the intercluster hopping 
matrix, respectively, and subsequent Fourier transfor- 
mation [3]. 

A substantial dependence of the system proper- 
ties on the band filling is expected for intermediate 
electron-phonon interaction, for which a large (ex- 
tended) polaron state is known to exist in the single- 
electron case. By contrast, for strong coupling, the 
polaron size collapses to a single site, so that different 
carriers will not overlap (interact). Hence, we take 
A = 1, the critical coupling for the large-to-small po- 
laron cross-over in the adiabatic regime uio/t < 1. Note 
that the abovementioned density effects are absent for 

LUo/t > 1 [1]. 

Figure 1 shows results for the one-electron spec- 
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Fig. 1. (color online) One-electron spectral functions A (k, u>) 
band fillings n. Here u> /t = 0.4, N = 10 and A = 1. (a) n = 0.1; 

tral functions A ± (k,uj) = — ImG ± (k,w) /tt, where G~ 
(G + ) corresponds to the (inverse) photoemission part 
of G, for fixed cluster size N = 10. 

In the low density case n = 0.1 (Fig. 1(a)), we 
can easily identify a (coherent) polaron band crossing 
the Fermi level Ef, the latter being situated at the 
point where A~ and A + intersect. The band flattens at 
large k, as known from single-polaron studies. Below 
this band, there exist equally spaced phonon excita- 
tions with small spectral weight, which reflect the Pois- 
son distribution of phonons in the one-electron ground 
state. Finally, above Ef, there is a broad incoherent 
band whose maximum follows closely the dispersion 
relation of free particles. 

As the density n increases (Figs. 1(a) -(c)) , the gap 
between the polaron band and the incoherent excita- 
tions is reduced, and the low-energy phonon peaks be- 
gin to broaden, until they have ultimately merged into 
a broad band for n = 0.4 (Fig. 1(d)). More impor- 
tantly, for n > 0.3, a polaron band can no longer be 
identified as incoherent excitations lie very close to the 
Fermi level. Finally, at n = 0.4, the spectrum consists 
of a broad main band of effective width ~ 6i cross- 
ing Ef, as in a metallic system, and low-energy ex- 
citations are available. This indicates that the pola- 
ronic quasiparticles, existing at small n, have dissoci- 
ated into dressed electrons due to mutual interaction. 
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(red dashed lines) and A^(k,w) (black solid lines) for different 
(b) n = 0.2: (c) n = 0.3; (d) n = 0.4. 

Therefore, the system can no longer be described in 
terms of small-polaron theory, as expected for interme- 
diate electron-phonon coupling at large n. 

To summarize, the density-driven cross-over from a 
polaronic system to one with moderately dressed elec- 
trons has been studied by means of CPT, whose unlim- 
ited momentum resolution significantly simplifies the 
interpretation compared to previous calculations [1]. 
Finally, the improvement of the cluster diagonalization 
represents an important step toward future studies of 
even more complex problems. 
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